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1. Introduction 

Here we discuss TMQCD at maximal twist with Nf = 2 degenerate sea quarks using config- 
urations from ETMC [|l|-|3|]. In particular we focus on neutral mesons, which have some unusual 
properties in TMQCD. Thus we evaluate the disconnected contributions which are needed for a 
study of neutral mesons We present a new method which, for twisted mass, allows many 
disconnected contributions to be evaluated very efficiently. 

We then present results for pions (exploring the charge splitting between tt" and 7r+). We also 
present results for flavour singlet pseudoscalar mesons (the 772 meson), for vector mesons (charge 
splitting and decays), and for iiavour-singlet scalar mesons. 

2. Disconnected diagrams - variance reduction 

TMQCD has a degenerate pair of m, d quarks with fermion matrix: Mj, j = it ijiy^ where 
Mw is the Wilson-Dirac matrix. Hence 1/M„ - l/Mj = -2//i(l/Mrf)75(l/M„). 

Consider the disconnected contribution /M„ — 1 /M^) where X is some 7-matrix and/or 
colour-matrix and the sum is over space. The conventional method involves solving 0^ = (1 /Mi,)^r 
with stochastic volume sources ^r, then ^X/M„ = L(i^*^0)r where the average is over noise 
samples (labelled r). However, the case mentioned above can be evaluated efficiently using the 
'one-end-trick' [^, ^. Then the required disconnected loop is given by 

This has signal/noise V / Vv^ = 1 which is much more favourable than the conventional 
method with signal/noise 1 / ^/V (here V = L^T). 

For example, at j3 = 3.9, = 0.004, L = 24, (where M{7i) ^ 300 MeV, a = 0.086 fm. La = 2. 1 
fm) taking X = i which is appropriate for the T]2 correlator, we evaluate the momentum-zero loop 
at a time-slice. We find a standard deviation of a = 1 8 from inherent gauge-time variation whereas 
the stochastic noise is a = 87 from 24 samples of volume source (conventional method) but only 
a = 7.5 from above method (with 12 samples). So 12 inversions give the disconnected correlator 
from all t to all t' with no significant increase in errors from the stochastic evaluation. 

For cases where this method cannot be used, eg. with VA75T3 i/a — > x^X maximal twist, we 
use hopping parameter variance reduction ^ instead. 

3. Pion order(a^) effects 

For ETMC data with M{k+) w 300 MeV, a = 0.086 fm and Lc? = 2. 1 fm, we show the ratios of 
correlators in fig. |l[ The disconnected pieces are seen to be relatively large - and reduce the charge 
splitting as found previously [Q]. The tt" is lighter than the 71+ , unlike a previous preliminary study 
of dynamical fermions [^. From a 4 x 4 fit to these correlations, we obtain Ama = 0.027(7). Since 
we expect rQ(m(7r*')^ — m(;r+)^) = c{a/rQY, we compare this expression with results from several 
lattice data sets [Q] in fig. ^ 

We see that, as expected, the flavour splitting decreases as a^. The sign and behaviour are 
consistent with Chiral PT and the nature of the phase transition [jl|] where m(7r") =0. We can use 
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Figure 1: Ratio of correlators between neutral and charged pions. The left plot is for local operators xIX 
for tt" and ^75;i; for 7r+ and the right plot for local operators XjAjs'^-iX for and XYaX for The curves 
give the ratio arising from the mass difference determined by the full fit. 



this determination to estimate the consequences of smaller lattice spacing, for instance less than 
20% pion splitting for m{Ti^) = 200 MeV provided ro/a > 8.2. 



4. Flavour singlet PS meson: 772 

In QCD, the flavour singlet pseudoscalar meson acquires a mass tlirough the anomaly, so 
is not a goldstone boson. It is important to check that this feature, which is linked to topological 
charge fluctuations, is reproduced in lattice evaluations. With A/y^ = 2 degenerate quarks, the flavour 
singlet pseudoscalar meson (called 172) is related to the experimental T]'(958) and is expected 
to have a mass around 800 MeV. We fit the T72 correlators (2x2 matrix with local and non-local 
operator \lfy5'^f ^ X'^iX at maximal twist) for f-range 3-10 with 2 states. We compare results from 
TMQCD [§] with older results (see ref. [|T^ for a review) in fig. |. Note that the ETMC results 



are at substantially smaller quark masses. We see that the T]2 mass is consistent with a constant 
behaviour in the chiral limit with m(T]2) ~ .88 GeV {rQm{r]2) = 2). 

We now discuss why the errors are so large for the T]2, despite the fact that we measure all t and 
t' , we use many gauge configurations and stochastic errors are small. The origin of the problem 
is that the signal for the disconnected part of the correlator comes from only a small part of the 
total data sample. For instance (at /i = 0.004 with 48 f-values for 888 gauge configurations) with 
I f' — f I = 10, 2.1% of the data contributes 26% of the signal. Thus the statistical impact of the data 
set is smaller than expected since parts of the data have big fluctuations (in a fermionic loop related 
to topological charge density). So even more configurations are needed to get reliable and small 
errors in the case of disconnected contributions. 
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Figure 2: Pion charge splitting 



5. Vector mesons 

We compare the local-local correlators (including disconnected parts for the neutral meson) for 
vector mesons in fig. Note that the disconnected contribution is negligible for the vector coupling 
to neutral p -mesons. We find that the ratio is consistent with constant (ie no mass splitting) and the 
value of that constant can be related to renormalisation constants. For the vector coupling to neutral 
p-mesons this implies (Z^/Zi/)^ « 1.5 at the finer lattice spacing, consistent with [|l2|]. We find 
agreement with a ratio of 1.0 for the tensor couplings, as expected since there is only one tensor 
renormalisation which then cancels. 

From fits to a 2 x 2 matrix of correlators (from the connected neutral contribution only here in 
f-range 8-18) we obtain p° masses (in lattice units). We also report values for the p+ masses [||] 
(from fits to a 4 X 4 matrix of correlators). These values are consistent with no flavour splitting for 



vector mesons as expected [11]. 
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am(p'') .400(25) .395(17) .419(17) .372(29) .346(12) 
am{p+) .416(12) .404(22) .428(8) .337(20) .337(12) 

We now consider decay transitions from the vector meson to two pions following the methods 
usedinref. [[131]. The transition p"(0) ^ 7r+(l)7r^(—l), where ^(tt) = ±2n/L, even for the lightest 
pion (M(7r+) ^ 300 MeV, L = 24a, a = 0.086 fm) is not open (note, however, that decay for a 
meson with non-zero momentum is open). For p''(0), there is an energy splitting of Ama = 0.19 
assuming that the two pion state has twice the energy of a 71+ with appropriate momentum. 

We show the normalised result in fig. H where it is compared to a two-state model [13|. The 
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value of the transition amplitude x is consistent with the empirical p decay width. One can also 
estimate the effect of this mixing transition on the p mass using the two-state model. This gives 
a downward shift of ma of .02 (eg from .41 to .39 for the p mass). This shift, induced by the 
proximity of the lightest two pion level, is comparable to our statistical error in determining the 
p mass. This suggests that we do not yet see major modifications of our p meson masses from 
mixing with the decay channel. 



6. Flavour singlet scalar mesons 



There is considerable confusion in allocating the experimental flavour-singlet scalar meson 
(/o) spectrum to specific content: since scalar glueball, uu + dd, ss, Tin and/or KK in an S-wave, 
etc., can all contribute. Lattice QCD can help considerably here, but it will be difficult as we now 
illustrate. 

We took a first look in TMQCD (here M{k+) w 300 MeV, a = 0.086 fm. La = 2.1 fm) and 
made a 2 state fit (? -range 6 to 23) to the 6x6 correlator matrix (P, 5, A4 both local and fuzzed 
at sink and source, including disconnected contributions). We find states ma = .103(5) (tt") and 
m'a = .227(28) (/o, energy consistent with 2m{n^)). 

Thus in the scalar channel we find a clear signal - but at the mass of two pions. This is not 
unexpected - but emphasises the problems of studying scalar mesons with light quarks in dynamical 
lattice gauge theory, where the light two-body state will dominate the correlators. 
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Figure 3: 772 mass versus quark mass. Here we summarise results with light pions and a < 0.1 fm (ro/a > 
4.5). 
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Figure 4: Ratios of correlator for neutral p-mesons to charged. The dotted lines guide the eye in the case 
that there is no mass splitting. 



7. Summary 

TMQCD allows efficient evaluation of disconnected contributions using a powerful variance 
reduction method. 

The n charge splitting goes to zero like as expected, and the sign is consistent with the 
nature of the phase transition. 

The flavour singlet pseudoscalar meson (772) has been studied to lighter quarks than previously 
and is consistent with a mass of around 800 MeV in chiral limit. 

For vector mesons we find negligible charge splitting and the decay (transition to 7i n) is 
accessible to study. 

Our study of the flavour singlet scalar meson (/o) with light dynamical quarks, finds the ex- 
pected n K contribution which will obscure further study of heavier states. 
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Figure 5: Transition k'^k' from the lattice (with M(7r+) « 300 MeV, a = 0.086 fm, La = 2.1 fm ). 

The solid line is from a two state model (p at rest and two pions with momentum q — ±2n/L) with energy 
gap Ama — 2aE{n) ~ am{p) ~ 0.19 and transition amplitude x/a. 
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